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Abstract. The energy bands and the global density of states are computed for superconductor / normal- 
metal superlattices in the clean limit. Dispersion relations are derived for the general case of insulating 
interfaces, including the mismatch of Fermi velocities and effective band masses. We focus on the influence of 
finite interface transparency and compare our results with those for transparent superlattices and trilayers. 
Analogously to the rapid variation on the atomic scale of the energy dispersion with layer thicknesses 
in transparent superlattices, we find strong oscillations of the almost flat energy bands (transmission 
resonances) in the case of finite transparency. In small-period transparent superlattices the BCS coherence 
peak disappears and a similar subgap peak is formed due to the Andreev process. With decreasing interface 
transparency the characteristic double peak structure in the global density of states develops towards a 
gapless BCS-like result in the tunnel limit. This effect can be used as a reliable STM probe for interface 
transparency. 

PACS. 74.45.+C Proximity effects; Andreev effect; SN and SNS junctions 



1 Introduction 

The artificial S/N superlattices consisting of alternating 
superconductor (S) and normal-metal or semiconductor 
(N) layers have been studied for some time already 
see also |10lllj . The recent advancement of nanofabrica- 
tion technology and experimental techniques |12| . as well 
as intrinsically layered structure of high-T c superconduc- 
tors |13ll4llftj has reinvigorated the long standing interest 
in understanding the effects inherent to clean supercon- 
ducting heterostructures |16I17I18| . The size and coher- 
ence effects have been studied recently for double barrier 
SNS and NSN junctions in the clean limit based on the 
solutions of Gor'kov and Bogoliubov-de Gennes (BdG) 
equations |19I2()I21I22I28I24| . 

In this paper we extend the previous approach of Tan- 
aka and Tsukada jS] and Plehn et al. [Sj to the more 
general case of superlattices with finite interlayer trans- 
parency. We present comprehensive and systematic anal- 
ysis of the influence of interface transparency on the quasi- 
particle band structure and density of states for wide range 
of the superlattice parameters. Due to the phase coher- 
ence of electronic wave functions the energy spectrum 
is gapless in superlattices with thin S layers and trans- 
parent interfaces [21], and splits into almost flat bands 
(transmission resonances) with decreasing transparency. 
For thick S layers, the subgap bands are formed due to 
the Andreev reflection [25] which leads to the conversion 
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of Cooper pairs in superconducting layers into correlated 
electrons and holes in the normal layers. Whereas the cal- 
culations are performed in the clean limit, the influence of 
impurities on the density of states can be taken into ac- 
by replacing the superconducting coherence length 
with an effective one, as shown by Halterman and Vails 
PH] in comparison with experiments of Moussy et al. ^2] ■ 
Our results for density of states in superlattices with layer 
thicknesses smaller than the superconducting correlation 
length, qualitatively confirm main features previously ob- 
tained by Bulaevskii and Zyskin 0] and Buzdin et al. [5] 
for atomic-scale layered systems within the tight binding 
approximation. 

2 The model 

The system under consideration is an S/N superlattice 
in the clean limit, consisting of alternating superconduct- 
ing and normal-metal (or semiconductor) layers of thick- 
ness ds and d/v, with insulating interfaces modelled as 
thin potential-energy barriers. The superconducting lay- 
ers are characterized by constant pair potential Zio, and 
zero phase difference, (f> = 0, is assumed across the super- 
lattice. Effective band masses and electrostatic potentials 
of the two metals are ms (tojv) and Us (Un), respectively. 
The superlattice is uniform in the x — y plane and the z 
axis is perpendicular to the layers. 

Quasiparticle propagation in the superlattice is de- 
scribed by the Bogoliubov-de Gennes equation 



ffo(r) A(t) 
A*(r) -HS(r) 



!?(r) = m{r), 



(1) 
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where &(r) = (u(r), t>(r)) T is the two-component wave 
function in the electron-hole space, the quasiparticle en- 
ergy E is measured with respect to the chemical potential 
/x, and the hamiltonian Hq within the superlattice period 
a = d]y + ds, for z £ (—(In, ds), is given by 



2m(r) 



+ S(z + d N )} 

+ C/(r)- M . (2) 



The first term is the quasiparticle kinetic energy in the ef- 
fective mass approximation |16I27| . the second term, with 

W = h 2 kFsZ/2ms, describes finite transparency of S-N 
interfaces modelled as (5-function potential barriers, and 
dimcnsionlcss parameter Z measures the barrier strength. 
Fermi energies in N and S layers are Efm = h 2 k 2 FN /2niN — 
Ll — Un and Eps = fi 2 k FS /2ms = fi—Us, respectively. We 
define the corresponding effective chemical potentials as 

= Miv(fc||) = e fn - {h 2 k 2 /2m N ) and ^5 = (is(k\\) = 
Efs — {h 2 k 2 /2ms), where k|| is the conserved quasiparti- 
cle momentum parallel to the layers. 

The pair potential A(r) should be treated self-consis- 
tently. For the sake of simplicity we used the stepwise 
model with A(r) equal to constant Aq in S and zero 
in N layers 8 19 . However, for S/N superlattices with 
thin S layers, the effective Aq can be taken as the space- 
averaged self-consistently determined pair potential, cor- 
respondingly smaller than the bulk value |23| . For super- 
lattices with thick S films, Aq can be set to the bulk value. 

Solutions of BdG equation, 



*(r) 



u(z) 
v(z) 



(3) 



in N and S layers can be written in the form 



«co) = Cl sin( ^ z) (0) + ° 2 cos{k " z) (0) 



and 



u{z) 
v{z) 



Here, fl 



C 3 sin(k N z 



= C5 sin(kgz) 



+ C-; sin(fc s z) 



C4 cos(k N z) 



(4) 



Ve 2 



A 2 k 



0' ™N 



l)+Cscos(k s z)(l 



V2m N (fi N ± E)/h 2 , k 



(5) 



v / 2ms(/xs ± f2)/H 2 and the BCS coherence amplitudes are 
u = V(l + 0/E)/2 and v = V(l - Q/E)/2. 

Complex coefficients C\ through Cs are determined 
from the boundary conditions at interfaces z = and 
z = —d^ inside the primitive cell 



u N (0) 
v N {0) 



u s (0) 
v s {0) 



1 fn' N m,k FS fu(0)\_ 1 fu' s (0)\ () 
^ U^v(O) ) + ~ S A U(0) ) ~ ^T s U(0) J ' {1) 



Ka ( u N{~d N ) 

v N (-d N ) 



iKa 



u' N (-d N ) 



itin \ v n(~^n) J m s 



us(d s . 
vs(ds] 



a f U N (-d N ) 

\y N (-d N ) 

= J_ hWs) 
m s \v' s (d s ) 



(8) 



(9) 



Here, the Bloch condition &{x,y,z + a) = e lKa ^(x,y,z) 
is used and the crystal momentum K is taken within the 
first Brillouin zone, K £ (— 7r/a, n/a). 

Dispersion relation E = E n ^.. (K) can be written in 
the following implicit form |Hj 



cos(Ka) = - A/4 ± yf (I?i/4) 2 - L> 2 /4 + 1/2 

= F ± (E,k ll ), (10) 

where D\ and D2 are defined in terms of dimensionless 
quantities E/A , k\\/k F s, Z, djv/Co, <W£o, m N /m s , 
Efn/Efsi and A /Eps (see the Appendix). 

Global density of states (for both spin orientations) 
per unit area of the cross section L x L y , averaged over a 
primitive cell, is given by 



g(E) = -±- S(E-E(k h K)) 

y (7,fe|| ,k 

= ~ I dk \\ k \\ E Y I dK *HE-E(k h K*)) 
n J i=+,- n J 

= ^J dk \\ fc n E 



i=+, 



dK 1 



dE{k h Ki) 



E(k«,K i )=E 



where if ± (E) are the solutions of Eq. (|1(J|) , and 





1 


darccos^E.fcu)] 


dE{k h K l ) 


a 


dE 



(11) 



(12) 



In accordance with Eq. (|10|l . the integration over k\\ [or 
I^S = MsC^ll)] i n Eq. i|llf) is limited to the intervals given 

by 

(L>i/4) 2 -i) 2 /4+l/2>0 (13) 



and 



|i^(£,fc||)| <1. 



(14) 



In the following, g(E) is normalized to the normal-state 
' value g = (msdskFS + ^w^TvfcFAr) /n 2 h 2 . 
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Fig. 1. Energy bands as a function of the S layer thickness 
ds, for S/N superlattices with thick N layers, djv = 3£o, trans- 
parent interfaces, Z — 0, and kii = 0. Andreev bound states 
[E < Ao) and geometrical resonances (E > Ao) for the cor- 
responding SNS and NSN trilayers are shown for comparison 
(solid curves). 
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Fig. 2. Characteristic dispersion of energy bands illustrated 
for djv = 3£o, Z = 0, ky = 0, and (a) ds = 3.0005djv, and 
(b) ds — 3.0013djv. Bands displayed in (a) are double degen- 
erate for all K. Dotted curves represent dispersion for finite 
transparency Z = 0.5. 

3 Energy bands and density of states 

The dispersion relation, Eq. i|l(J|) . is solved numerically 
and the global density of states is calculated from Eq. i|TT)l 
for various superlattices and for zero phase difference (j> = 
0. In the following, we focus on the influence of finite in- 
terface transparency on quasiparticle band structure and 
density of states. For simplicity, this is illustrated for equal 
effective masses and Fermi wave- vectors, mjv/ms = I 
and kpN/kFS — 1- Superconductors are characterized 
by the bulk value of the pair potential Aq/Efs = 1CP 3 , 
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Fig. 3. Energy bands as a function of Z, for ky = and for 
two particular S/N superlattices: (a) ds = djv = 3£o and (b) 
ds = 15£o, dN = 3£o- Arrows indicate the bound states in the 
tunnel limit. 

which corresponds to the zero-temperature BCS coherence 
length £ = ti 2 k FS / (irm s A Q ) - 10 3 A. 

Energy bands for S/N superlattices with thick N lay- 
ers, transparent S-N interfaces and quasiparticles prop- 
agating perpendicular to the layers (k|| = 0) are shown 
in Fig. H Quasicontinuum of energy states corresponding 
to the crystal momentum within the first Brillouin zone, 
K €E [—it/a, Tr/a), is indicated by shading the band width 
calculated from Eqs. ifT^jl and 

For the corresponding SNS trilayer, Andreev bound 
states, E < Aq, in the normal interlayer of thickness d/v, 
for zero phase difference across the junction, transparent 
interfaces, and k|| = are given by j^H] 

i\ , 1 A\l 1 , 1t .\ 

— — = 7r raH — arccos — — - — — , 15) 
A L 7T \A J\d N /£o y ' 

where n = 0, 1, In this case the Andreev bound states 

are double degenerate. Geometrical resonances, E > Ao, 
for the corresponding NSN junction with S interlayer of 
thickness ds, and Z = 0, kii = are given by 

7T = \l l + n2 7T7TY2> ( 16 ) 

Ao y {ds/ior 

which follows from the condition ds(kg — kg) — 2nn, 
where n — ±1,±2, .... At these energies the Andreev re- 
flection vanishes and the electron is transmitted without 
creation or annihilation of Cooper pairs 22 23 24 . Both 
Andreev bound states and geometrical resonances of the 
corresponding SNS and NSN trilayers are shown in Fig.^ 
for comparison. 
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In S/N superlattices with thick S layers, the energy 
band structure above the gap, E > Aq, is also affected by 
the Andreev process |17I18| . With increasing d,N, the band 
structure dependance on ds/diq remains qualitatively the 
same as in Fig. with compression and lowering of en- 
ergy bands that enter the superconducting gap (221 - An- 
dreev reflection is the fundamental mechanism that de- 
termines the quasiparticle band structure in S/N super- 
lattices. However, qualitatively the same results as shown 
in Fig. are obtained for semiconductor / normal-metal 
superlattices |cSUI31| . Characteristic dispersion of energy 
bands, shown in Fig.0 is illustrated in Fig.0for two close 
thicknesses of the S layer. For some layer thicknesses the 
energy bands are double degenerate for all K, Fig. [3(a), 
in contrast with the usual degeneracy at high-symmetry 
points only (at the center and the ends of the first Brillouin 
zone) , Fig. (b) . These two types of dispersion alternate 
rapidly with the change of layer thicknesses on the atomic 
scale kp 1 , while the band width changes on the macro- 
scopic scale. 

Finite interface transparency, as well as mismatch of 
effective masses and Fermi wave- vectors [22] , lift the de- 
generacy in E(K), Fig. [3 and change the band structure, 
Fig. [21 For large Z, energy bands split into pairs of flat 
bands independent of K, and there is a significant change 
of the band energy below the superconducting gap. Ap- 
proaching the tunnel limit for Z» 1, pairs of adjacent flat 
energy bands transform into bound states of isolated films 
defined by dsk s = n^n and d^k^ — n-i^. However, this 
does not imply that the energy band splitting and decrease 
of the band widths due to the flattening will be visible in 
the E vs. dg(jv) plot (cf. Figs. and EJ. Energy levels 
for Z ^S> I oscillate rapidly with layer thicknesses on the 
atomic scale kp , Fig. [SJ so that E vs. d S ( N ^ curves fill the 
energy space quasicontinuously on the macroscopic scale, 
Fig. This implies erasing of the band structure and lo- 
calization of quasiparticle states in real superlattices with 
finite interface transparency and slightly unequal layers. 

Previous analysis has been made for quasiparticles that 
propagate perpendicular to the layers. Dependence of en- 
ergy bands on k\\ , i.e. on the effective chemical potential 
£is(fc||), is illustrated for Z = in Figs. El (a) and 0(a). 
Band widths decrease with the increase of , and bands 
split into pairs of bound states (flat bands) for very large 
parallel momentum [H], similar to the tunnel limit. The 
effect of erasing the band structure with finite interface 
transparency is enhanced with the increase of k\\ . Corre- 
sponding changes of the global density of states are shown 
in Figs. El (b) and (b). Integration of Eq. (|TT|) is per- 
formed over the shaded regions in Figs. El (a) and0 (a), 
where Eqs. lfH?|) and (fHfl are satisfied. 

Now we shall focus on energy bands and the density of 
states in thin-layer S/N superlattices, where coherence ef- 
fects are pronounced and ballistic transport is more likely 
to take place [SIJ. Dependence of energy bands on the 
superlattice period is illustrated in Fig. for d$ = djv, 
k|| = 0, and for both Z = and Z = 1. It can be seen that 
the band structure in transparent thin-layer superlattices 
differs significantly from the thick-layer case considered 
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Fig. 4. Energy 'bands' for djv = 3£o, ky = and Z = 4. 
Shading is produced by the rapid oscillatory dependence of 
flat energy bands on the S layer thickness. 




Fig. 5. Rapid oscillatory dependance of energy bands on the S 
layer thickness for nontransparent superlattices. Energy disper- 
sion E(K) is shown by shading for all K, solid curves represent 
E(Q). 



in Ref. 0. For thin layers, dispersion of energy bands is 
significant, with only a small part of the lower band lay- 
ing below Aq. Energy bands as a function of feii, and for 
various interface transparencies are shown in Fig. El For 
Z = and d$ — g?at, the most striking feature is the onset 
of the lowest energy band at the midgap, practically for 
any k» . This is not the case for thick-layer superlattices, 
where band energy decreases more rapidly down to zero 
with the increase of feii, resulting in the left-side "tail" of 
the subgap peak in the density of states, Figs. El and 
For thin-layer superlattices, finite interface transparency 
introduces the resonance effect: energy bands penetrate 
periodically below the midgap with the increase of k». 
This is more pronounced as Z gets larger, Fig. El The cor- 
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Fig. 6. (a) Normalized energy bands E/Ao as a function of 
MS = Mst^u), for S/N superlattice with ds = 15£o, rfjv = 3£o 
and Z = 0, and (b) the corresponding global density of states 
ff(-E) normalized to the normal-state value g. Global density 
of states for Z — 0.5 (dashed curve), and in the tunnel limit 
(dotted curve) are given for comparison. 
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Fig. 7. (a) Normalized energy bands E/Ao as a function 
of us ~ (J,s(ku), for S/N superlattice with ds ~ djv = 3£o 
and Z — 0, and (b) the corresponding global density of states 
g(E) normalized to the normal-state value g. Global density 
of states for Z = 0.5 (dashed curve), and in the tunnel limit 
(dotted curve) are given for comparison. 



responding global densities of states for various interface 
transparencies are shown in Fig. 1101 

For transparent interfaces, the density of states is BCS- 
like with the energy gap E g smaller than the pair potential 
Aq. The value of E g for transparent interfaces, equal ef- 
fective masses, and equal Fermi energies can be obtained 
from the well known dispersion relation 



cos[(K ± ± k Z F)a] = cos(q6 ds) cos(q (In) 

— S^ 1 sm(qS ds) sm(qdjy), (17) 



which is a special case of Eq. (EH. Here, k zF = VfefT 



5 = H/E, and q — mE/h 2 k Z F- For ds, d N 
Jill and (JUJ exactly follows [T] 



/• 2 
II ' 



E„ 



J,. 



1 + d N /d s ' 



0, from Eqs 



(18) 



Practically, this simple relation remains valid for the layer 
thicknesses up to one or two coherence lengths £o: due 
to the weak variation of the bottom of the lowest energy 
band with the layer thickness, Fig. 0(a). 

With decreasing interfacial transparency, the subgap 
peak in g(E) at E g decays, and the usual BCS coherence 
peak at A reenters as the superconducting layers become 
more isolated. In the tunnel limit, the BCS peak at Aq is 
completely restored, Fig. ^| (dotted curve in the bottom 
panel). For thicker layers ds ~ d^ ~ £qj the coherence 
effects are less pronounced and the tunnel limit behav- 
ior is practically reached for smaller Z ~ 1. Previously, 
this double peak structure in the density of states of S /N 
superlattices is obtained within the tight binding approx- 
imation for atomic-scale layered systems, and apparently 
observed in high-T c intrinsically layered superconductors 

mm- 



4 Conclusion 

We have derived the dispersion relation for superconduc- 
tor / normal-metal (semiconductor) superlattices in the 
clean limit, generalizing the previous expression of Plehn 
et al. [S] to include an arbitrary interface transparency and 
mismatch of effective band masses. The obtained general 
dispersion relation is used for numerical analysis of the 



< 




Fig. 8. (a) Energy bands for S/N superlattices with ds — dN, 
k|| = and Z = 0. (b) Erasing of the band structure with 
decrease of interface transparency is shown for Z — 1. 
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Fig. 9. Energy bands as a function of us ~ fis{ku), for S/N 
superlattice with thin layers, ds = djv = 0.1£o, and for various 
interface transparencies. 



influence of interface transparency on energy band struc- 
ture and density of states in metallic S/N superlattices. 
Although we used stepwise approximation for the pair po- 
tential, our results will not be altered significantly by the 
fully self-consistent numerical calculations |S] , if an effec- 
tive pair potential (smaller than the bulk value) is taken 
for thin S layers, and simply the bulk value in the cases of 
thick S layers, low transparency, and mismatch of Fermi 
wave-vectors or band masses. Our results confirm previ- 
ously obtained features in the metallic S/N superlattices 
|6I7I8| . including the limiting cases of double barrier SNS 
or NSN trilayers |M1^0li>ll^l^liM| . and are in a good 
qualitative agreement with the results obtained within the 
tight binding approximation |4I5| . 

Consequences of the quantum interference effect are 
strong and rapid (on the atomic scale) geometrical os- 
cillations with layers thickness of the energy dispersion 
in transparent superlattices, and of the almost flat en- 
ergy bands (transmission resonances) in the case of finite 
transparency. Oscillations in the latter case are practically 



w 






■ Z = 50 







E/A„ 



Fig. 10. Global density of states for S/N superlattice with 
thin layers, ds = djv = 0.1£o, and for various interface trans- 
parencies. Tunnel limit is indicated in the bottom panel (dot- 
ted curve). Note that the effective Aq varies with Z, being the 
smallest for Z = (top panel) and reaching the bulk value in 
the tunnel limit (bottom panel). 



within the band width of the corresponding fully trans- 
parent superlattice. Practically, this could imply the lo- 
calization of quasiparticle states in superlattices with low 
interface transparency. 

Characteristic changes of quasiparticle band structure 
with decreasing interface transparency are suitably re- 
flected in the global (averaged on the lattice period) den- 
sity of states, which can be directly measured by STM 
techniques |14ll2j . Oscillations of the density of states 
are simply related to the band structure for transparent 
superlattices with thick layers 6 . However, superlattices 
with the period smaller than the coherence length, when S 
and N layers lose their individual properties due to strong 
phase coupling by Andreev scattering, do not differ sig- 
nificantly from a bulk BCS superconductor, except for a 
subgap peak in the global density of states instead of the 
superconductor coherence peak at Aq For transpar- 
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ent interfaces, position of the subgap peak is simply re- 
lated to the lattice parameters. For finite interface trans- 
parency, we find the characteristic double peak structure 
in the global density of states With decreasing trans- 
parency the subgap peak decreases, only slightly changing 
the position, while the coherence peak at the effective Aq 
grows, and the density of states develops towards the gap- 
less BCS result for the bulk superconductor in the tunnel 
limit. We point out that this double-peak structure of the 
global density of states in small-period clean-metal S/N 
superlattices can be used as a reliable experimental probe 
for interface transparency. 

We are grateful to Ivana Petkovic, Milos Bozovic, and Boris 
Grbic for useful discussions. The work has been supported by 
the Serbian Ministry of Science, Project No. 1899. 
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Appendix 

From the boundary conditions, Eqs. the dispersion relation, Eq. (fTT)|) . is expressed through L>i = Di/D and 

F> 2 = D 2 /D 0l where 



Dq =m^ k N k^ k s kg (u 2 — v 2 ) 2 , (19) 

D 1 =F + F 1 (Zk F8 ) + F 2 (Zk FS ) 2 1 (20) 

D 2 =G + Gi{Zk FS ) + G 2 (Zk FS ) 2 + G 3 (Zk FS ) 3 + G 4 {Zk FS ) 4 . (21) 

Here, m r — mjv/ws, Fq through F 2 and Go through G4 are given by 

Fq =m r (u 2 — v 2 ) I [k^ kg Sg (fc^ 2 + ml kg 2 ) + k^ kg s^ Sg (fc^ 2 + m 2 r kg 2 ) ] u 

[fcjV fcg Sjy Sg (fcjv ~t~ kg ) + fcjy kg Sjy- Sg (fcj\r + 7Yl r kg ) j V 

+ 2m r k^k^ kg fcj [v 2 (c^Cg + c^cj) ~ ^ (44 + c n c s) ]}■. ( 22 ) 

F 1 = - 2 m 2 . (it 2 - u 2 ) jfc^ fc^ [kg Sg (c^u 2 - c^v 2 ) + kg sj (cj w 2 - c^u 2 ) ] 

+ m r kg kj [kjf s% (4 ^ - c s « 2 ) + k N $n {c~s ^ - c% v 2 ) ] } , (23) 



and 



F2 = — rr? r (it 2 — v 2 ) [k N s^ (fcg s^ u 2 — kg s s v 2 ) + k~^ s N (kg s s u 2 — k s Sg w 2 ) ] , (24) 

Go = <Sjy Sjy Sg Sg (fcjv ^iV ^5 ) ^ ) 

- 2 (u 2 - u 2 ) |m r k^ kj, [k% c^ s^ {ks c s 4 " 2 ~ 4 4 s s « 2 ) + k N 4 (4 4 s s ^ - kg Cg Sg~ v 2 ) ] 

+ ml kg kg [k N c N s^ (fcg c s Sg u 2 — k s Cg s s v 2 ) + kjj c^ s N (fcg Cg s s u 2 — kg c s Sg v 2 ) ] | 
+ m 2 (2 k N k^ {sg sj u 2 v 2 (c^ c^ - l) (fcg 2 + fcg 2 ) 

+ fe s 4 [ (l + 2c^c^Cg 4) (" 4 + ^ 4 ) - 2uv 2 (cg4 +cj i cj i + c^ct 1 cg4)]} 
+ s jv s jv {2&s kgiiv (c s — l) (fcjv + fe^ ) 

+ s s Sg[ii (fcjv kg + fc^ fc s )+« (fc^ fc s + fc^ feg ) ] }1, (25) 

Gi = — 2 m r (it 2 — u 2 ) 

\ + ( -f+7+ - T7-2-2 .+2-2 . 2 , — 2 /-2 -2\1 

x Sjv I m r Sjy jc s fcg s s [Kjy it — k N v + m,. fe s [u — v J J 

-,- + T/+2-2 ,-2-2 . 2 ,+2 f -2 -2\ 1 l 

+ c s kg Sg [kj, u — k N v + m r kg (u - w ) ) 
+ c jv ^iv { ~ 2 Cg cj m 2 fcg fcg (it 2 — v 2 ) + Sg Sg [k~^ 2 (it 2 — v 2 ) + m 2 . (kg 2 u 2 ~ kg 2 v 2 ) ] }J 

-(- Cjy fejv ^ 2 Cg 171 r /bg (eg 711 r Sjy /bg -(- Cjy kjy Sg) (it U ) 

+ Ss{-2c N 4 m r fc JV fc S (^ 2 _ + S N S S [ k N 2 (" 2 — " 2 ) + m r ( fc S ^ ^ ~ k S ^ ^) ] } ) ] J ( 2 6) 



G2 = — m 2 (it 2 — u 2 ) 



^2cg m r kg {2cg m r s^ s^ fcg (— u + v 2 ) + Sg [ — c^ k^ s^ (u 2 -2v 2 ) + c^ k^ s^ (-2 m 2 + u 2 ) ] } 
+ Sg ( - 2 c 4 } fc^ [eg" m r s^r fcg (ii 2 - 2 v 2 ) + 2 c^ fc^ Sg (u 2 - ii 2 ) ] 

+ s jv { s jv [ k N 2 + k N 2 + m r { k s 2 + 4 2 ) ] s s (^ 2 _ " 2 ) + 2c^ 4 m '- fc iv 4 (~2ii 2 + v 2 ) }) j , (27) 



G 3 = 2mf. (u 2 - n 2 ) 2 {m r fcg s^ s 4 ^ s s Cg + Sg [fcjy Cjy s 4 ^ s s + s N (m r k s c s + c% s s ) ] }, (28) 



Mihajlo Vanevic, Zoran Radovic: Quasiparticle states in superconducting superlattices 



9 



G4 = rn\ s N s s Sg (y 2 — v' 2 ) , (29) 

where = sm(k^dN), = cos(k^dN), s;| = sin(fcg ds), and Cg = cos(fcgds). For Z — and m r = 1 expressions 
for D\ and Z?2 reduce to the results given in Ref. [H]- 



